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1 Introduction 



The phenomenon of vacuum polarization is a triumph of flat space Quan- 
tum Electrodynamics (QED). The enhanced high energy coupling strength 
it predicts has not only been verified experimentally, it also provided the 
crucial paradigm for understanding renormalization group flows. In spite of 
these successes, the enhancement is only about 2% at the highest energy 
accelerators we have so far been able to build pp. 

Loop effects such as vacuum polarization derive ultimately from the re- 
sponse to virtual particles (in this case, electrons and positrons) which are 
the quantum field theoretic manifestation of 0-point motion. It has long been 
realized that these effects must be strengthened in an expanding universe, 
essentially because the expansion of spacetime tends to hold virtual quanta 
apart [2]. Leonard Parker made the first quantitative computations in the 
late 60's [3]. He found that the effect is largest for massless particles that 
are not conformally invariant jl] , which includes massless, minimally coupled 
scalars and gravitons [5]. At a fixed expansion rate the effect increases with 
the cosmological acceleration. During primordial inflation massless virtual 
quanta are actually ripped out of the vacuum, which is thought to be the 
origin of primordial scalar and tensor perturbations [6]. 

From this discussion one can infer that the largest possible vacuum polar- 
ization occurs during de Sitter inflation (which has the highest acceleration 
consistent with stability) and derives from massless charged particles that 
are not conformally invariant. These conditions are realistic. The measured 
value of the scalar power spectrum, and the current upper bound on the 
tensor-to-scalar ratio [TJ [5] suggest that the deceleration parameter (which is 
minus the acceleration) of primordial inflation was less than -.993 [9], which 
amply justifies taking the de Sitter limit of -1. Primordial inflation is thought 
to have occurred at such an enormous scale that all known charged parti- 
cles would have been effectively massless. However, massless fermions are 
conformally invariant on the classical level, which means they can only ex- 
perience the expansion of spacetime through the conformal anomaly. This is 
responsible for the vacuum polarization of ordinary QED being only slightly 
enhanced during primordial inflation [10J. 

Much larger effects can come from a charged scalar (such as the compo- 
nents of the Higgs which become the longitudinal parts of the W in low 
energy physics) provided it is massless and not conformally invariant. Davis, 
Dimopoulos, Prokopec and Tornkvist made the remarkable proposal that 
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such a particle might even endow the photon with mass during primordial 
inflation [TTJ [T2] . Their idea was confirmed with a dimensionally regulated 
and fully renormalized one loop computation of the vacuum polarization 
from Scalar Quantum Electrodynamics (SQED) on de Sitter background [T3] . 
Although the one loop effect grows without bound [14], a nonperturbative 
resummation of the leading secular terms reveals that the photon mass ap- 
proaches a constant value of about 1.8163 times the Hubble constant [T5] . 

The quantum-corrected, linearized Maxwell equations in an arbitrary 
metric g^ u read, 



0, 



gg up g^F pa {x)] + Jd A xfw](x; x')A u (x') = J»{x) . (1) 



Here A p (x) stands for the vector potential, [^TF^a;; x') for the vacuum po- 
larization and J M (x) is the current density. It is immediately obvious that 
the same vacuum polarization that reveals the effective mass of dynamical 
photons can also be used to compute the electrodynamic response to point 
charges and current dipoles. That is the purpose of this paper. 

Throughout this paper we employ a spacelike metric. Section 2 derives 
corrections to the field strengths of point charges and dipoles from the one 
loop vacuum polarization of flat space SQED. That serves as a useful cor- 
respondence limit and also illustrates the basics of the more challenging de 
Sitter computation. In section 3 we review the de Sitter geometry. Section 4 
presents the actual de Sitter computation, with the messy details consigned 
to two appendices. Our conclusions comprise section 5. 



2 Vacuum Polarization in Flat Space 

The purpose of this section is to work out the effect at one loop from a mass- 
less, charged scalar in flat space. Our analysis parallels recent computations, 
made using the Schwinger-Keldysh formalism [T6l [T7] , of the one loop effect 
of gravitons on electromagnetism [18] and of a massless, minimally coupled 
scalar on linearized gravity in flat space background [19]. We begin with 
some considerations following from the general tensor structure of the vac- 
uum polarization and from the fact that we only have the result to some 
finite order in the loop expansion. Then the one loop correction is derived 
for a point charge and for a point magnetic dipole. 
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2.1 General considerations 



Poincare invariance and transversality constrain the vacuum polarization at 
any order to take the form, 



"IT 



(x;x') = &-dTf v -d»d v Il(x-x') 



&-&jf v -&^& v Tl(x-x') . (2) 



This means we can integrate by parts to express the effective field equation 
© as, 

d u F v ^(x) - Jd 4 x' U{x-x')d' u F^ l {x') + (Surface Terms) = J»{x) . (3) 

The causality of the Schwinger-Keldysh formalism precludes there being any 
spatial surface terms, or any surface terms at the upper limit of the temporal 
integration [20J. There can be surface terms at the lower temporal limit, 
which corresponds to the time when the state was released [21]. We shall 
assume that these are completely absorbed by perturbative corrections to the 
initial state wave functional [22] ■ If we agree to exclude these state corrections 
from n(x — x'), the linearized effective field equations become just, 



d u F u »{x) - Jd 4 x'U(x-x')d' u F^{x') = J"{x) . 
The one loop result for II [x— x') is [T5] , 



n (1) (x-x') 



ad 4 



96tt 2 



#(At-Ax)jln[/i 2 (At 2 -A£ 2 



(4) 



(5) 



Here d 2 = T]^ v d^d v = — <9 2 + V 2 , At = t — t', Ax = \\x — x'\\, and a = e 2 /A-n is 
the loop-counting parameter of SQED. Because we do not possess the higher 
order corrections, equation (jlj) can only be solved perturbatively. That is, 
one expands the field strength in powers of a, 



P(0) , p(l) _i_ ei(2) , 



/ll/ 



fJ,V 



(6) 



and then distills (jlj) into terms of the same order. Because the current density 
is zeroth order we have, 



d v Fj${x) 
d v F${x) 



d 4 x'U^(x-x')d' u Fi,°J{x' 



(7) 
(8) 
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and so on. We can therefore define the source of the one loop field strength 
to be the "one loop current density," 



jjP(x)= jd A x'U.^{x-x')J^x') 



(9) 



It remains to explain that we will always solve for the field strength tensor 
directly, rather than going through the intermediate step of finding the vector 
potential. Note that contracting the Levi-Civita density into the gradient of 
d a F av gives, 



c pa„UQ gap 



Substituting ([9]) into ([8]) and then using (fTO]) implies, 



(10) 



(ii) 



This form is advantageous in view of expression ([5]) for IT^^x — x'). Up to 
a homogeneous solution, we can write the one loop field strength tensor as 
FIm = dn<?v ~ ®v3u, 1 > where jjp is obtained by simply removing a factor 
of d 2 from 



ad 2 



96vr 5 



d A x'0(At-Ax)l\n fi 2 (At 2 -Ax 2 ) -1 }J„(x') . (12) 



The ambiguity regarding homogeneous solutions can be settled by appealing 
to initial conditions. 



2.2 Response to a point charge 

The current density of a stationary charge q at the origin is, 

J"(f, x) = q5^5 3 (x) ^ J^t, x) = -qSy(x) . (13) 
Taking the initial time to be tj and inserting f JT3|) in expression f fl2|) gives, 

ft— X I 



J^(t,x) 



d 2 db' < In n(t-t'-x) +ln /i(t-t'+x 



1 



96vr 2 



967T" ,/tj 

<9 2 <! -2xln(2 / ux)+3(x-t+t / ) + (t-t/-x) lnL(t-tj-a:) 



(14) 



+(t—ti+x) In fi(t — ti + x 



(15) 



Before acting the d'Alembertian it is useful to specialize it to functions of t 
and x, 



d 2 



X 



dix 



^(d t -d x )(d t +d x )(t-x)-^(d t -^ x )(d t ^d x )(t+x) . (16) 
2x 2x 



We now use ( JIB]) to act the d'Alembertian in ([T5 



907T Z X X 



t — ti — X 



(17) 



The homogeneous terms would vanish if acted upon by another d'Almbertian, 
but they also drop out when tj is taken to the infinite past, 



lim jU(t,x) 

tj — ^ — oo ^ 

The one loop field strengths follow, 



aq5° ln(2/xx) 



24vr 2 



x 



aq x 
24vr 2 x 1 



ln(2/ix) — 1 



6n 



ln(2/ix)-ll x F$\t,x 



F^(t,x) 



. 



(18) 

(19) 
(20) 



From expression f fT9|) we see that the regular izat ion scale /i gives rise to a 
characteristic length, L = e 1 /2fi. Outside this length the charge g is screened 
by quantum corrections, whereas it is enhanced for x < L. Why this happens 
becomes clearer from examining the one loop current density, 



Jfo(t,x)=d 2 J(i ) (t,x) 



24tt 2 \" 



+ oo x 5 3 (x) 



x" 



(21) 



Of course it is the positive contribution at the origin which enhances the 
electrostatic force for x < L, and the negative cloud of charge density screens 
it for x > L. 

Vacuum polarization can be understood physically from the response of 
virtual charged particles (in this case massless scalars) to the classical source. 
The energy-time uncertainty principle allows a pair of virtual particles with 
mass m and wave number k to exist for a time At ~ 1/ \/m? + k 2 . During 
this time the partner whose charge is opposite to q will be pulled into the 
source, and the partner with the same charge as q will be pushed away. By 
itself, that would lead to a negative induced charge at the origin, which is 
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exactly opposite to ( l2~Tj) . However, two features complicate the physics of 
vacuum polarization in this model: renormalization and the masslessness of 
our scalars. The first feature means that what we call the one loop current 
density J^(t,x) actually includes an infinite constant times the classical 
current density. Had the mass been nonzero, this infinite constant would 
have been chosen to completely null the one loop correction to the fields at 
spatial infinity. However, the masslessness of our charged scalars means that 
there continues to be an effect out to arbitrarily large distances. In that 
case, one chooses the renormalization constants to make the one loop field 
strengths vanish at some fixed distance, L = e 1 /2/x. 



2.3 Response to a point magnetic dipole 

Shrinking a current loop with magnetic dipole moment m down to the origin 
gives the following current density, 



x 



J u = 



J 



-mx W 3 (f) . (22) 



This differs from the current density f fT3|) of a point charge only by the 
replacement, 

q5$ — ► -e 0tip °m p d a . (23) 

Because the derivative must be spatial it can be partially integrated, then 
converted from V' to —V, so the result for J~^(t,x) follows by making the 
same replacement (I2"3"j) in expression f llSp . 



lim JL (t, x) 



a 



tj^-oo - W ' ' 247T 2 

It is useful to 3 + 1 decompose (1241) . 
J^\t,x) = 0, 



^ pa m p d a 



ln(2//x) 



x 



jj 1] {t,x) = ln(2/ix)-l 



a 



6tt 



The resulting field strengths are, 



x 



4nx 3 



F 0i {t,x) 

Fij (t , 3?) 







ae^ k ( (3x k x-fh — m k 



24tt 2 { 



X" 



ln(2/ix) — 1 



+ 



(m —x x-rh) 



x° 



(24) 

(25) 
(26) 

(27) 
(28) 
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In view of the relation = —e l i k B k we can read off the one loop correc- 
tion to the magnetic field from 



BM(t,x) = -^[ln(2H-ll x flC)(t,i) - - x l"'-**'"*' j (29 ) 



6tt 

where (away from x = 0) the classical result is, 



Qtt Attx 3 



m(t,x) = {3 ™'™~™ ] . (so) 



Expression (1291) represents the same sort of screening we found in (1191 for a 
point charge. The second term in ( )29|) is just the residue necessary to keep 
B transverse with the coordinate-dependent screening factor. Except for the 
usual divergence at the origin, the one loop current density which supports 
this field is, 

j) ^^ = -^-2— (^0). (31) 
The current rotates clockwise with respect to m. 



3 de Sitter Geometry 

The de Sitter geometry is the unique, maximally symmetric solution of Ein- 
stein's equation with a positive cosmological constant. Any coordinatiza- 
tion x^ of de Sitter can be described by a mapping x^ — > X A (x) to the 
4-dimensional submanifold of 5-dimensional Minkowski space such that, 

- (X ) 2 + (X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 = r lAB X A X B = -L • (32) 

Here H is the Hubble constant, which is related to the cosmological constant 
by A = 3H 2 . The de Sitter metric in x^ coordinates is given by, 

dX A dX B 

9M = ^T^ab • (33) 

The most convenient coordinates for our work are open conformal coor- 
dinates, x^ = (r],x). Whereas each of the spatial coordinates x l runs from 
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— oo to +00, the conformal time runs from 7/ — > —00 (the infinite past) to 
77 — > (the infinite future). The 5- dimensional embedding is, 

X° = + (34) 

X i = ax 1 , (35) 



X 



1 



" 1-H 2 (x 2 - V 2 ) , (36) 



2H 

where the scale factor is a = —l/Hrj. The inverse transformation is 



____ JvabXW X* 
1 #(X°+X 4 ) ' " H{X°+X 4 ) ' 1 J 

From expression ( 1331) we see that the invariant element is, 

ds 2 = a 2 (r]) —drj 2 + dx-dx , a{r]) = ——— . (38) 
l -I ' Hf] 

Hence = a 2 r]^ u . 

Because 4-dimensional electrodynamics is conformally invariant we can 
reduce many aspects of the de Sitter computation (in conformal coordinates) 
to familiar manipulations in flat space. One feature of conformal coordinates 
that sometimes disturbs the mathematically inclined is that they do not cover 
the full de Sitter manifold. One can see this by adding f )34|) and fl36|) to get, 

X° + X 4 = ^ > . (39) 
H 

This is of no import because rj = const defines a Cauchy surface, so infor- 
mation from any part of the full de Sitter manifold must enter the conformal 
coordinate submanifold as an initial condition. The full de Sitter manifold is 
irrelevant if one regards de Sitter as a paradigm for the geometry of primor- 
dial inflation. From that perspective the conformal coordinate submanifold 
is just a special case of the larger class of spatially flat, Friedman-Robertson- 
Walker (FRW) geometries which are relevant to inflationary cosmology. 

We shall also be interested in static coordinates, x^ = (t,r). Although 
the time coordinates runs from t — > —00 to t — > +00, the spatial radius obeys 
< r < 1/H. The 5-dimensional embedding is, 



x° 


= isinl 






X 1 


= r\ 


X 4 


= Icos 







■H 2 r 2 , (40) 

(41) 

-H 2 r 2 . (42) 
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The inverse is, 



t = itanh" 1 (^) , r'= , £— = ■ (43) 



HJ VAB XAXB 



Combining with f l3"3"|) we find the invariant element 

dr 2 



ds 1 



1-H 2 r 2 



^ + i ^72 2 + r W • ( 44 ) 

1 — H z r z 



Of course this is why the system is called "static coordinates." 

As we have mentioned, mathematically oriented physicists sometimes 
imagine that subtle errors must result from failing to formulate quantum 
field theory on the full de Sitter manifold. Static coordinates pose an equal 
peril for confusing those who seek to deny the possibility of secular effects as- 
sociated with the continuous creation of massless, minimally coupled scalars 
and gravitons during inflation. Quantum field theory cannot be formulated 
in static coordinates because t = const does not constitute a Cauchy surface. 
The static coordinate t is the proper time of an observer in free fall at the 
origin. Although no point outside the static coordinate system can causally 
influence this observer, that is only true for r — 0. All other points in the 
static coordinate system are in causal contact with points from outside. 

The preceding considerations preclude computing the vacuum polariza- 
tion in static coordinates, but they in no way prevent us from transforming 
results from conformal coordinates into static coordinates. If the source is 
at the static coordinate origin then an observer at fixed r would appear, in 
conformal coordinates, to be moving towards the origin with precisely the 
velocity needed to counteract the expansion of spacetime and maintain a 
constant physical distance from the source. Such an observer will experience 
effects associated with boosting the conformal field strength, and other ef- 
fects associated with shrinking the conformal coordinate separation to zero. 
By substituting (I4"0ll4"21) in ( 1371) we see that the coordinate transformation is, 

_ g "~ ~~ ia ^ ^ 

71 ~ Hy/l-H 2 r 2 ' X ~ Vl-H 2 r 2 ' ^ 
The various components of the Jacobian are, 

dr/ e~ Ht dx l Hr z e~ m 



dt y/l-H 2 r 2 ' dt Vl~H 2 r 2 ' 
dr/ Hrie~ Ht dx l 8 l ie~ m // 2 ;•'/■'( ' 



(46) 



+ z ... ,3 • (47) 



dri [l-# 2 r 2 ]i ' dr* ^\-H 2 r 2 [1-H 2 r 2 \2 
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We denote the static coordinate field strengths with a tilde, 

e -2Ht ( \ 

e -m 2 r 

Fa = \-Hr l F 0j + H^F 0t + (l-H 2 r 2 )F tJ 



(48) 



+H 2 r i r k F kj -H 2 r j r k F ki \. (49) 



The factors of e seem to drive F pu to zero but one has to keep in mind 
that a factor of 1/x 2 grows at fixed r, 



Jim 



x 2 r 2 



1-H 2 r 2 



Note also that the scale factor is, 



a = e H Wl-H 2 i 



(50) 



(51) 



4 Vacuum Polarization in de Sitter Space 

The plan of this section is the same as that of section 2. We begin with 
some general considerations in which the one loop result for the [ At n i/ ](a;; x') 
is presented in open conformal coordinates, and the one loop currents it 
engenders are evaluated as much as possible for a general source. We then 
specialize the classical source to a point charge and a point magnetic dipole. 
In each case the one loop current density and the one loop field strengths are 
first derived in open conformal coordinates, and then transformed to static 
coordinates. 



4.1 General Considerations 

An important simplification associated with conformal coordinates in D = 4 
spacetime dimensions is that all the scale factors in the classical Maxwell 
equation cancel out, 



-gg» p g^F pa 



= d v 



(52) 
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This means much of the flat space formalism of section 2 still applies. In 
particular, the equation for the one loop correction to the field strength is, 

d 2 F$ = d,jV - ft, j<D , (53) 

where d 2 = r]^ u d^d u = —d 2 + V 2 is the flat space d'Alembertian and the one 
loop current is, 



Jj»(x) = T) W d 4 x' \m» {1) ](x-,x')AW(x') . 



(54) 

The retarded solution for flB^l) is, 

F$)(x) = -ftx> 6(/ f~ A ^ X) ^4V)-^V)] , (55) 

where the conformal coordinate separations are Arj = r\ — rf and Ax 
r x'\\. 

Because the massless, minimally coupled scalar is not de Sitter invariant, 
the vacuum polarization it engenders contains two distinct tensor structures. 
One of these is proportional to the covariant transverse projection opera- 
tor, whereas the other one is proportional to the purely spatial transverse 
projection operator constructed from rf v = rf v + 5q(5q|j 



(x\ x') 



&„d a U F (x; x') + U c (x; x') + U G (x; x') 



+ \f]^f} pa -fJ flp r} va ]d' p d a U K (x-,x') . (56) 

It is convenient to break the covariant structure function up into the old, 
flat space contribution lip, a term He like the conformal anomaly, and a 
nonlocal de Sitter contribution He. At one loop order one has [US, E]> 



11$ (x; x') 
ng } (x;x') 
ng } (x;x') 



a 



96vr 5 
a 



d 4 



e(Ar)-Ax)\ln \x 2 (Arj 2 - Ax 2 



— — ln(a)5 4 (x — x') 

D7T 



aH 2 a 



d 2 



a'6{Ari-Ax){ ln\H' 2 {Ar] 2 -A 



x 



+ 1 



(57) 
(58) 
(59) 



^^One might wonder if there is any advantage to employing a covariant representation 
based on covariant derivatives and de Sitter invariant basis tensors. However, a systematic 
investigation of this formalism reveals that it is cumbersome and that it obscures the 
essential physics |23) . The procedure for converting our noncovariant — but simple - 
representation to the complicated and counter-intuitive covariant representation can 
be found in [Ml. 
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The one loop noncovariant structure function is 

e(AT]-Ax)\\n\H ,2 (A7] 2 -Ax 



_ ,Tj4: n 2 n /2 ( 
Yl^(x-X') = u, \ .. \ ..> ) 1,, U'2; \ ,2 \ ,2, 



+ 2 



(60) 



Here H' = e 7_1 if , where 7 ~ 0.57721 is Euler's constant. Although the spa- 
tial dependence of the various structure functions is limited to the coordinate 
separation, x — x', factors of a and a' complicate the temporal dependence. 

The various one loop current densities can be expressed in terms of the 
classical current density J M and the classical field strength Fff), 



ad 4 



96n 2 
aH 2 



dV^(Ar/-Ax)|ln[/i 2 (A^-A^ 
ln(o)F®(x)l = ^-\-hx{a)J^{x)+HaF ( 



l\Ja(x') 



6tt 



(o). 



X 



8tt 2 



ad 2 Jd 4 x' a'6(Ar]-Ax) 

H%A V 2 -Ax 2 )]+i\f^{x') 



x < In 



Jl K {x) 



aH a 
4ir 2 



4 o 2 r 

-^-/c/Va' 2 #(Ar/-A.x) 
i 2 J 



x < In 



H' 2 (Ar ] 2 -Ax 2 )\+2\d' j Fjl\x') . (64) 



(61) 
(62) 



(63) 



Of course Jh (x) is the same as the flat space result of section 2. Note 
that Jq K (x) = because the index ft is purely spatial. It is useful to 3 + 1 

aH 2 a 



decompose J} G (x) 



Jl G {x) 



Jl G (x) 



8tt 2 
aH 2 a 



d 2 d 4 x'a'9(Ari-Ax)\\n H' 2 {A V 2 -Ax 2 ) +l\J {x') , (65) 



or d x' a'9(Ar) — Ax) I In H' 2 (A v 2 - Ax 2 ) +l\J l (x') 



+ 



aH 3 a 2 



d 2 ld 4 x'a'9(Ar]-Ax){ In H' 2 (Ar] 2 -Ax 2 ) +1 \F^\x' 



aH 3 a 



s _ 2 d 2 Jd 4 x'a' 2 e(Ar ] -Ax)^\n[H' 2 {Ar ] 2 -Ax 2 )\+l { jF^\x') . (66) 
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We should also comment on the range of coordinates. Unlike the flat 
space case, it is not possible to release this system in free vacuum infinitely 
far back in the past because perturbation theory breaks down after a finite 
interval of conformal time [15J. The one loop results ( 1571160]) were computed 
with initial time rjj = —1/H, corresponding to aj = 1. For a source at the 
origin, the factors of ^(A^ — Ax) in ()57H60p imply that there is no one loop 
effect for radii x outside the range, 

rji < 7]—x =>• Hx < 1 — . (67) 

a 

For future reference we define the scale factors a± evaluated at conformal 
times r] ± x, 



1 _ 1 _ at ~L±Hr 
a± = ~H{n±x) ~ a~ l THx~ e Vl?#r ' (68) 

The late time behavior of these scale factors depends critically upon whether 
or not one fixes the conformal radius x or static radius r, 



M — > =Ft|- , (a ± ) ^e H \^. (69) 

V ^Jx fixed Hx ^ /r fixed \j \ 

4.2 Response to a point charge 

The electromagnetic coupling to a point charge q with worldline x M ( r ) ls i 

S q = qJdrrA,( X ). (70) 

Because S q is independent of the metric, the current density of a point charge 
is unchanged from (|13p . nor are the 0th order field strengths changed from 
section 2.1, 

J»^x) = q5»5\x) , 4°W) = 7-4 > F$\ v ,x)=0. (71) 

Air x A 

One consequence is that J^ F must agree with expression (1211) . and the asso- 
ciated one loop field strengths must agree with ( 11911201) . 



Jl F M = {x ± 0) , F^( V ,x) = -£[ln(2/ix)-l] x F<$(r, 



x 



(72) 
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Substituting (1711) into (1521) gives, 

fc)^M^) , ^ M = ^ = -*{^}. 

(73) 

Hence the C-type charge density represents a fractional screening of the 
classical charge density by — |Mn(a), with the current density J}° carrying 
off the positive charge to infinity We can easily solve for the C-type electric 
field by noting, 

5 2 F f ( V ,x) = d j} c ( V ,x)-d^ c ( V ,x) = atJ^MzL). (74) 



' ; \ 247r 2 ' 



x 



Therefore, the C-type field strengths are, 

F^(v, 2) = ~ x F<$ ( V , x) . (75) 

Substituting (1711) into (1641) implies that all components of the If-type 
current density vanish. It is better to approach the G-type current density 
and its associated field strength indirectly. Suppose the temporal and spatial 
components of the current density take the form, 

J 1G (T7,f) = f(rj,x) , Jl G (r ] ,x) = d l g(v,x) . (76) 

Current conservation implies, 

= -d v f(r],x) + V 2 g(rj,x) g(r] 1 x) = J da/(l-—)x?d rl f(r] 1 a/). (77) 

The magnetic field is obviously zero and we can find the electric field by 
noting, 

d 2 F^( V ,x) = d v Jl G { V ,g)-diJ^ G ( V ,^ , (78) 
= ftj^/-/}, (79) 

= -<9 2 <9 4 [ X dx'(l--)x'f{ri,x') . (80) 
Jo V X ' 

Hence everything follows from the zero component, 

Jl G ( v ,x) = ^ £dx'x' 2 d v J^ G (r],x') , (81) 

Ft G (r,,x) = -^£dx'x' 2 J l G ( V ,x'), (82) 
Fl G ( V ,x) = 0. (83) 
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Expression (j82p is recognizable as the integral form of Gauss's law, and ex- 
pression (1811) is the Maxwell displacement current. 
It remains to evaluate Jq G (^,x), 



jf( V ,x) 



^f^V P drj'a'\\n\H'(AT]-x)\ +\n\H'(Arj+x)] +1 1. (84) 

87T Jr\i I I 



Expression (1841) is sufficiently intricate that the best strategy is to treat each 
of the three integrands separately, and to make an additional distinction 
between the two derivative operators which result when the d'Alembertian 
is specialized to functions of just r\ and x as in ([IS]) . 

d 2 — ► ^(^-^X^+a^-x) - L^-d^dn+d^ri+x) . (85) 

The first of the integrands in (1541) depends only upon 77 — x so it vanishes 
when acted upon by the first operator of (155]) . 



^(^-a^+d^-x) f Jr/a' In [iT(A77-x) 



. 



(86) 



A nonzero result emerges when this same first integrand is acted on by the 
second operator in (1551) . 



~(d v -d x )(d v +d x )( V +x) f V dr/a'lnlH'i&n-x) 

ZX Jr\i L 



1 



X 

2a- 

x 



Tj — X 



— {d v -d x ) / drj'a'ln H'(Arj-x) , (87) 



m 



| 7 -l-ln(a_) + ln(l-^)|. (81 



(We remind the reader of the scale factors a± defined in (jBSJ.) Acting the 
first part of the d'Alembertian (155]) on the second integrand in ( 154]) gives, 



^-(d v -d x )(d v +d x )( V -x) f dr)'a'ln\H'(Ar)+x) 

2x Jr)i 1 



—(d ri -d x )(r}-x){-ln(a-)-ln(2Hx)+ln(l — - 
x v cu 



2a, 



x 



ln(a.)-ln(2^) + ln(l-^) + I-^- 



(89) 
(90) 
(91) 
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The second part of the d'Alembertian (I85p acts on this second integrand to 
produce, 



— (d^-d^+d^rj+x) I d V 'a'ln H'(A V +x) 



Tj — X 

VI 



- — (d 11 +d x )(r l +x)\n(2Hx), (92) 

x 

2a. 



x { 



\n(2Hx) + 



2 2aHx 



(93) 



And the final term of the integrand in 
bertian can be kept together, 



is simple enough that the d'Alem- 



rj—x 



d 2 / d V 'a' = d 2 



VI 



ln(a_ 
H 



2a. 



x 



(94) 



Summing expressions (156"]) . (}9"Tj) . (I9"3"j) and (19~4"|) . and multiplying by 
the prefactor of (154"j) gives the G-type charge density, 



J l G (rj,x) 



aqH 2 a ( a. 



4tt 2 



x 



a + 
+— 

x 



ln(a_)-ln(2ifx)-m 1-- 
-ln(a_)-ln(2#x)+m 



1 

1 

a. 



1 



2 2a#x 

1 1 

2 ~ 2aifx 



(95) 



Expression ( 195]) is complicated but its import becomes clear at late times, 
either at fixed x or at fixed r = x/a, 



1G 




x fixed 



{ 2ln (^)- ln ( 1 -^ V ) +0 ( 1 )}' 



1G\ 
o J 

x=r/a 



aqHa 

4lT 2 X 2 

aqH 2 a 3 f 21n(a) 2#rTn(l+#r) 21n(2i2r) 



47r 2 r l + #r 



(96) 

+ 0(1)). (97) 



In each case we see that the one loop charge density (which is minus Jq G ) 
is opposite to q, implying screening. We also see that the effect from Jq G 
is much stronger than that of Jq° (powers of a versus In (a)), which is itself 
stronger than that of Jq F . 

The G-type contribution to the one loop electric field follows from substi- 
tuting each of the ten terms of expression (1951) in (|82|) . The integrals are not 
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illuminating and have been consigned to expressions (115011159]) of Appendix 
A. After considerable rearrangement the final result is, 



a 



7T 



x F^\r], x) x I aHx 



(l-7)[l-ln(- + #a; 
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aHx 



In 



1\2 



-In 



1 + aHx 



21n(a)--ln 



1 + aHx 



+ l- 7 -ln(2) 



-Li 2 
2 



4-Li 2 



1 1 
1 

2 2 



-aHx 



12 2 



ln 2 (2) + Li 2 









'—Ex 


— Li 2 


-+Ex 






-a J 



-ifV 

l-a 2 #V 
(9 



The symbol Li 2 (z) denotes the dilogarithm function, 



U 2 (z) 



ln(l-t) ^ z k 1. , 7T 2 ^ I 

(X6 



^fc 2 2 1r2( 2) 6 § /, 



fc=i 



(99) 



k=l 



Note that we have employed the identity, 



vr 2 1. 



Li 2 (-x) -Li 2 (l-z) = -ln(x)ln(l + x) - y - ^Li 2 (l- r 



(100) 



Expression (jUSj) is unwieldy. One can understand it better by taking the 
limit of late time at fixed x, 



x fixed 



x x aEx x J 7 (Ex) 

ix An x 6 



(101) 



The proportionality function is positive definite, 



?{Hx) = l-MHx) 4-2 In 

L J Ex 



1 + Ex 



1-Ex 



ln(l-ifV) . (102) 



(Recall that 7 ps 0.57721 is Euler's constant, so 1 - 7 0.42278.) For small 
if x it has the expansion, 

T-(ifx) = (I-7) 4- (1-7) + \(Exf + ^(Exf + ... (103) 



The largest Ex can get is one, at which point, 

jr(i) = 3-7-2 ln(2) « 1.03649 



(104) 
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Combining expressions ( T72|) . (175]) and (110 ip gives the total one loop cor- 
rection to the electric field at late times for fixed co-moving position x, 



7 ar 

0t , 



x fixed 



--xF ',°». 

7T 



6 



ln(2/is)-l + -ln(a) + a#a;xJ r (Fz) \. (105) 
J 6 



It will be seen that the 1G contribution totally dominates the 1C and IF 
contributions. The one loop correction we have found is consistent with 
screening from the vacuum polarization induced by the vast ensemble of 
charged scalars produced during inflation. The presence of a scale factor 
means that the one loop correction cancels the tree order effect after only 
about ln(a) ~ ln(7r/a) ~ 6 e-foldings! Of course this does not mean fixed 
x observers see the field of a negative charge after that time; what happens 
instead is that they see effectively zero field strength. Perturbation theory 
breaks down because screening has become nonperturbatively strong. 

The result is curiously different at fixed r in static coordinates. To see 
why, let us first evaluate the 1G contribution ([98"]) at late times holding 
r = aHx fixed, 



r fixed 



a q a r , . . 

■- x — x m(aj 

7r Att r A 



-2la(l+Hr) + (l-j)Hr . (106) 



This is negative definite for all < Hr < 1, which is consistent with anti- 
screening. Of course that is what happens in flat space as well, when one 
gets very close to the source. Note that the factor of ln(a) cancels between 
the IF and 1C contributions at fixed r, 



6 



ln(2yux) — 1 H — ln(a 
J 6 



ln(2//r)-l 



(107) 



Hence the 1G contribution is still dominant. 

To derive the electric field strength in static coordinates we must still 
transform the vector indices. Substituting (I106j) into expression (HB]) gives, 



(4 1 



a 



q r 



r fixed 



x x Ht 

(0) x ^ x Ht 



-21n(l+Fr) + (I-7) Hr ,(108) 



0; 



21n(l+ifr) - (l-7)#r 



(109) 



We see that the one loop correction at fixed static coordinate r (1 1 9 j) has 
both the opposite sign and a much slower growth than the result in con- 
formal coordinates ( 11051) . The physical interpretation seems to be that the 
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static coordinate observer experiences the logarithmic running of the electro- 
magnetic coupling which is built into the Bunch-Davies vacuum, no matter 
what coordinate system one employs. The effect becomes nonperturbatively 
strong after about Ht ~ 7r/a ~ 430 e- foldings. 



4.3 Response to a point magnetic dipole 

In conformal coordinates the current density and classical field strengths of 
a point dipole are unchanged from flat space, 



J°(r],x) = 
F^( V ,x) = , F^( V ,x) 



J(r],x) = —m x VS 3 (x) 



e ijk 

4tc 



m fc V 2 -m-Va 



1 

x 



110) 
111) 



It follows that the IF contributions to the one loop induced current density 
and field strengths are the same as we found in section 2.3, 



J l F (rj,x) 







j" F (r],x) 



a mxx 
3vr 2 x 4 



F \ F ( V ,x) = , Flf^x) 



IF i 



ae 



24tt 2 



m k V 2 -m-Vdi 



(112) 

-i^}. (113) 



Note that we have neglected delta function terms in (I112j) . 

Substituting (IllOtlllip in expression ( )62|) gives the 1C contribution to the 
one loop induced current density, 



J lc (r/,f) = 



Jl c ( V ,x) 



a ln(a) 

67T 



s ijk m j d k S s (x) 



(114) 



Using these currents, with a few partial integrations, in relation ( 15"5]) produces 
the associated field strengths, 



ijk 



24tt 2 



mPd k d 



jHHx+l) 



or 



24vr 2 



m k V 2 -rh-Vd k 



(\n(Hx+±) 
\ x 



115) 
116) 



Recall that we are assuming < Hx < 1 — 1/a. 

It is best to add the G and K currents together before inferring the field 
strengths they induce, so we derive them together. From expression flB3|) we 
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see that Jq G (i],x) vanishes because the classical charge density is zero for a 
point magnetic dipole. Substituting the other classical values ( IllOHllip into 
expression f )66|) and performing some partial integrations reduces the induced 
G-type current density to a single integral over conformal time, 



Jl G ( V ,x) 



aH 2 a 



in' 



r\—x 



8 2 e ijk m 1 d k / dn'a'lki H' 2 {A V 2 -Ax 



ni 



117) 



Of course the temporal component of J (t),x) always vanishes. From ex- 
pression ([M]) and the classical relation djF^\r],x) = —e^ k m^dkd 3 (x) we 
infer, 



aH^a 2 



An 2 



q—x 



VI 



e ijk m j d k / dri'a' 2 {ln H' 2 (A V 2 -x 2 ) + 2 . (118) 



The next step is to perform the conformal time integrations in expressions 
(11171) and (j!18p . For the G-type current density comparison with (184"]) reveals 
that we can read off the result by making the replacement q — > e K i m?d k in 



aH 2 a 
4tt 2 



e ijk m 3 d k { — 



x 



2Ex' 



)-ln 


r 1 ' 


1 


1 


1 


-7 1 

[ 2 






a_- 


2aHx 



+- 



X 



-ln(a_)-ln(2ifx) + ln 



1- 



1 i 1 



+ 



1 



2 2aHx 



;il9) 



For the K-type current density it is best to change variables from rf to 
a' = — 1 / Hrf in expression f 1 1 1 8 [) , 



J- (??,*) = 
aH 3 a 2 



aH 3 a 2 
' 4tt 2 ' 





r 1 1 - 




r 1 


1 i 


: J da' jln 




+ ln 


-a 7 






■a' a_- 







+ 2 7 , (120) 



47T 2 



e^ k mP d k \^{a + — a_) ln(a + — a_) — (a_ — 1) ln(a_ — 1) 



-(a + -l)ln(a + -l)-(l-a_)ln(a + ) + (l+a_)ln(a_)-2 7 a_ \ .(121) 



The final step before combining the G-type and f^-type currents is to 
put them in a common form by isolating distinct logarithms and making 
judicious use of the identities, 



2a+a_ = a(a + + a_) 



a_i_ — a_ , 



Hx 



aHxa± = ±a± =p a 



(122) 
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The results are, 
aH 2 



4tt 2 



4vr 2 



e ijk m j d k { (a+ +a_) ln(2#:r) - (a+ - a_) ln(a_) 

I X L 

-o_ln(l — -)+a+In(l — -)-(7+l)a_ 1 , (123) 
a+ + a_-2a) ln(2#;r) 



+ (a + + a_ + 2ai/x — 2a) ln(a_) + (a_ + aifx — a) ln( 1 



a i - 



aHx — a) ln(l )+27(a_ — a) 1. 

V a + J J 



(124) 



Adding expressions HI 23[) and (11241) gives a current density which is much 
simpler than either of them separately, 



4vr 2 



x 



-2a\n(2Hx) + 2(a^ + aHx-a) ln(a_) 



-a(l-Hx) ln(l ) +a(l+Hx) ln(l-— ) -a_+7(a_-2a) 



.(125) 



Finally, it is useful to introduce the symbol J gk (t], x) for H 2 times the part 
of ( 11251) within the curly brackets, 



J GK ( V ,x) = ^ k m j d k J GK ( V ,x) . 



An 2 



(126) 



From (I53"j) we see that the one loop field strength induced by (I125p obeys 
the equation, 



d 2 F GK ( V ,x) = ^Jf%f)-^jf%£) . 
The retarded solution ([55]) takes the form, 



where, 



F° K (V, 2) = d,J u GK ( V} x)-d v J GK ^ x) 



AnAx 



(127) 

(128) 
(129) 
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We can commute the differential operator in (I126p through the Greens func- 
tion to define a scalar function J GK {rj,x), 

J? K M = ^ k m^d k J GK (r,,x) , (130) 

—^ k m^d k fd 4 x' ^ (Ar/ ~ Ax) J GK ( V ', x') . (131) 
In 2 J AtcAx 



4vr 2 

In terms of J' GK (r],x) the electric and magnetic field strengths are, 

ijk 



F T(v,x) = -^mW k d J GK ( V ,x), (132) 



Air 



ae 



ijk 



F GK (r],x) = -^ r [m k V 2 -m-Vd k \J GK (r ] ,x) . (133) 

It remains to evaluate the function J~ GK (ri,x). Note that for a function 
f(r),x) of just i] and x = \\x\\, the integral against the retarded Green's 
function gives, 

- /dV ^~ AX) x f( v ', x>) = -i- fdrf [ X dx p x'f( V ', x') . (134) 
J 4ttAx 2x J Vi J\x-a v \ 

In our case the function f(t],x) is J (j],x) and we must remember that 
causality makes it vanish for x > rj — r/j . We can therefore express 
as the sum of three integrals, 

JG*( n,x) = - — II dff dx'+ df]' dx' 

v " I 2x {Jn+M^L 'J v - v '-x Jri-x Jx-v+v' 

+ Ljv'l dx'\x'J GK (r)',x').(135) 



x—rj+rj' 



These integrals are evaluated in Appendix B. Here we present just the leading 
results for large a with x fixed, and with r = x/a fixed, 

j£? dx ^ l ^-xG(Hx) , jgZ,^ ^ x K{Hr) , (136) 

X X 

where the functions Q(z) and JC(z) are, 

G(z) = 2 7 + 21n(2)-ln(^)-zln(l-^ 2 ) , (137) 
K(z) = ^+21n(l + ^)+Li 2 (^) . (138) 
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The function Q{z) is positive at z = so the GK terms dominate over 
the F and C terms at late times for fixed x. 



ae 



fixed x 



fixed x 



An 2 



m j d k d ' 



ln(a) 



ae 



4tt 2 



x g{Hx) 

ln(a) 



x 



x Q{Hx) 



(139) 
(140) 



Near x = the term in curly brackets is the positive constant 27 + 2 ln(2) 
times ln(a)/x. This screens the classical result of —1/x, although only loga- 
rithmically. 

For small z the function fC(z) has the expansion, 



/C(z) = 2\n(2)z-2z 2 + 0(z 3 



(141) 



The linear term drops out from the one loop field strengths, so the GK terms 
are negligible compared with the F and C terms for large a at fixed r. The 
resulting field strengths are, 



(i) 
0; 



ae 



'J k 



fixed r 



247T 2 



mPdhdt 



fc^0 < 



-\n{2nx)+\n(Hx+\) 



x 



ae 



ijk 



24vr 2 



m j f k x Ha 3 x 



H 2 



ae 



fixed r 

ae ljk 
24tt 2 



247T 2 



! (l + i/r) 2 
\n(2/j l x)+\n(Hx+l)\ 



x 



(142) 
(143) 
(144) 



x a x 



V 0,,r / 



2/ir 



m k — 3m-rr k 



+ ^7 



m k — 2fh-rr k 



m —3m-rr 



1+Hr 



1 ^ 


m 


—m-ff- 


r 


- (1 


+ Hr) 2 J 



(145) 



We must employ expressions ( HHl) and ( 1491) to convert (11431) and (11451) 
into the static field strengths. The final result for the (late time form of the) 
static electric field is relatively simple, 



7(1) 

Oi 



ae 



247T 2 



m 



j r k xHe H Wl-Hh 



ln( 



l+Hr 
2[ir 



H 



r(l+Hr) 



■ (146) 



The factor of e Ht makes this seem like a substantial modification but one 
should keep in mind that the tree order result carries the same factor in 



23 



static coordinates, 



7(0) 

Ot 



e ijk 

An 



w?f k X 



He Ht \fl — H 2 r 2 



(147) 



It is therefore better to view the one loop correction as a time independent 
screening of the classical result by the r-dependent enhancement (for small 
r) of about a/Qn x ln(l/2/w). This is exactly what one finds in flat space 
( |29|) and is a manifestation of the running of the electrodynamic coupling. 
The tensor structure of the one loop static magnetic field is complicated but 
its leading term for small r bears the same relation to the tree order result, 



',(0) 



',(1) 



Jjk 



x e H Wl-H 2 r 2 



An 



rm — 3m-rr 



at 



tjk 



2An 2 



x e 



Ht 



Vl-H 2 r 2 \\n 



1 + Hr 
2jir 



H 2 vm k — m-rr k 



1-H 2 r 2 



m 



■3rh-rr 



+ 



(148) 
(149) 



5 Discussion 



In this study we have explored the effect that the inflationary production 
of charged, massless and minimally coupled scalars have on the electric and 
magnetic fields induced by a point charge and by a point magnetic dipole. 
The sources were held fixed at the spatial origin. We derived results for the 
fields perceived by two observers, one stationary in conformal coordinates 
and the other stationary in static coordinates. 

The conformal observer sees the one loop correction (11051) to the elec- 
tric field of a point charge screen the classical result exponentially rapidly in 
physical time. That is to be expected because the physical distance between 
this observer and the source grows exponentially, and the intervening space 
is filled with a nearly constant density of charged particles. In contrast, the 
static observer perceives the one loop correction fll09[) to slowly enhance the 
classical result. It is more difficult to understand this result. A possible 
interpretation is that it derives from the running of the electrodynamic cou- 
pling which is built into the scalar vacuum. However, it should be noted 
that the enhancement does not come from merely probing the one loop cor- 
rection in flat space ffl9|) at the ever-smaller coordinate separations needed 
to keep the physical distance to the source constant. That effect is entirely 
canceled by the conformal anomaly term (1751) ; the enhancement originates 
in the intrinsically de Sitter contribution (|98p . 



24 



Our results for a current dipole are interestingly different. The confor- 
mal observer sees much weaker screening fll40p which grows only linearly in 
physical time. The static observer sees an enhancement (11491) which is con- 
stant in time but grows as one approaches the source. Because the static 
result derives from the flat space and conformal anomaly terms it is the most 
simply understood in terms of the flat space result (|2"9"|) . and is clearly a 
manifestation of the running of the electrodynamic coupling. 

As far as we know this work is unique in having solved the effective field 
equations for the constrained part of a force field on de Sitter background. 
There have been many previous studies of the propagation of dynamical 
quanta on de Sitter, both from the production of scalars [25J E2 ttH E51 EZ] 
and from the production of gravitons [2SJ EH ED]- There have also been 
many flat space background studies of how massless scalars affect the force 
of gravity [311 [19] , and of how gravitons affect electric and magnetic forces 
[32j ITS] . But no one has previously combined de Sitter and force laws to 
reveal the surprising dichotomy between screening in conformal coordinates 
to anti-screening in static coordinates. 

Although this seems to be the first time one loop corrections to force 
laws have been studied on de Sitter, it will not be the last. Computations 
have already been made of the one loop graviton self-energy from inflationary 
gravitons [33] and from inflationary scalar [31] , and the vacuum polarization 
from inflationary gravitons will soon be completed [35]. These results can 
be used to quantum- correct the effective Maxwell and Einstein equations, 
and then proceed to a study of the same type we have just completed. We 
anticipate that the current work will play an important role in organizing 
and understanding these studies. The project to quantum-correct Maxwell's 
equations to account for inflationary gravitons is especially noteworthy be- 
cause it is a straightforward approach to check the puzzling screening effect 
that Kitamoto and Kitazawa have inferred from graviton loop corrections to 
the off-shell effective field equations at fixed sub-horizon scales |36j. 
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6 Appendix A: Integrals from Subsection 14.2 



The following integrals (times — a /it x FqP) sum to give our result (198]) for 



F^ in section 4.2, 



aH 2 f X dx'x'a'_\n(a'_) = - ln 2 (a_)-^ ln 2 (a) + — ln(a_)-ln(a) + a#2 , (150) 
jo 2 2 a„ 

-aH 2 [ X dx'x'a'_hx(2Hx') = \n(l + aHx)\n(2Hx) + U 2 (-aHx) , 
Jo 



-aH 2 faVx'a'. In fl — - 
Jo v a_ 



—aHx\n(2Hx) + aHx , 
a — — ) lnfl — -Wa-f) lnfl-~ 

+Li 2 (^)-Li 2 (-^)+a#x 

/ 1 \ f x / 1 \ r a 

yy + -)aH 2 I dx'x'a'_ = (7+ _ ) ln(a) ln(a_) — aifx 



a// 2 / dx' 
Jo 



- ln(a) ln(a_ 

2aH 2 y 1 2 v 



Z" 33 / a \ a 
—aH 2 J dx'x'a' + \n(a'_) = hif^J ^(l - aHx)-\ ln(a_) — ln(a) + aifx , 

-Li 2 (^)+Li 2 (---afTx) , 

-aH 2 J X dx'x'a' + \n(2Hx') = - ln (^) ln(l-a#z)+a#zln(2#:r)-aifz 

+Li 2 (l)-Li 2 (l-aFx) , 

a# 2 f X dx'x'a' lnfl -4-1 = -fa-—) lnfl — -Wfa-l) lnfl-- 

Jo v a!. y v a + y v a + 7 v ' v a 



a + - 

-Li 2 (-)+Li 2 (— )+«fl.c 



-aH 2 f dx'x'a'. 



-aH 2 / dx 



o 2aH 2 



- \n(l — aHx) — -aHx 
- hi(l — aHx) . 



151) 

152) 
153) 
154) 

155) 
156) 

157) 
158) 
159) 



The integrands of (I155H1561) and (I158H159]) diverge at x' = —rj, which results 
in the ill-defined factors of ± ln(l — aHx). However, these divergences cancel 
between (115511156]) and between (I158H159]) . as do the problematic logarithms. 
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7 Appendix B: Integrals from Subsection 



We begin with some general comments. The expression (11351) we must evalu- 
ate is a double integration over three regions which can be termed "J" , "IJ" 
and "///", 



J GK (v,x) 



1 

2x 



<dr]'dx ; + / di]'dx'+ / dr( dx 1 
i Jii Jin 



x'J GK (r]',x') 



(160) 



Each of these double integrals takes the form, 



dr] dx , 

JB JD(r}') 



where the various limits are, 



A! = r]-x , A n = -(rji+rj+x) , A IH 



7/ 



If 

d = r(-T}i 
Dj = f]—x—r]' 



B n = Aj 
Cu = Ci 



Bni = Mi i 
Cm = rj+x-rj' , 



D 



li 



D 



in 



-Dj. 



(161) 

(162) 

(163) 

(164) 
(165) 



Of course the integrand is of crucial importance. From expressions (1 1 2 5 [) 
and (11261) we see that xJ GK (rj,x) can be broken up into seven components, 



Xi(?7,x 
X 2 (ri,x 
234(77, z 

X 9 (?7,x 



-2H 2 a 2 ln(2#x) , 
+2# 2 aa_ ln(a_) , 
+2H 2 a 2 (l-Hx) ln\-H(r)-x 

-a 2 (l-Hx)\n l + H(r]-x 

+a 2 (l + Hx) In l + H(r]+x 
(7 — l)aa_ , 
-2a 2 . 



(166) 
(167) 
(168) 

(169) 

(170) 
(171) 
(172) 



In each case the x' integration can be expressed in terms of elementary func- 
tions. For T10 the r/' integration also results in elementary functions. I2 and 
Ig give polylogarithms, about which more later. The remaining integrands 
- Zi, Z34, Z 56 and Zf 8 — all take the form of a' 2 times the x' derivative of 
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an elementary function of 77' and x' . For these integrands the best strategy 
is partially integrate on 7/ after having performed the x' integration, 

, l2 r c y\dF(i,x') 



H drj'a' 2 dx 
Jb Jdm) 



(r)') 8x' 



HjX'a'^F^', C(r)'j) -F(p', D(rf)) } , (173) 

= *' C) -F(rf, D) } [ -J "\'a'^ { F (rf, C) -F(t]', d) } . (174) 
It turns out that the surface terms in fll74p cancel between the three regions, 
a'^F^', C I )-F{r 1 , ,D I ) I =a Al ^F(A I ,7 1 -x- m )-F(A I , 0)^-0, (175) 

a' |f(t/, C n ) -F(rj', D n ) j = a AlI \f{a iu ^( v +x- Vi )] 
-F(A II ,^( Vl - V +3x))^-a Al ^F(A 



!,r]-x-r]i 









Am | 


a!< 


F{r ) \C III )-F{r 1 \D in ) 


\ 


= 0-a Al A 








Bui I 



■F(Aj,0)| ,(176) 
(177) 



-F(A n ±to-v+3x))\ . (178) 



Integrands Z 2 and X 9 , and the volume terms fll74p from I 1 , X 34 , I 56 and 
I78, all give rise to 7/ integrations of the form l/r/' times logarithms. Inte- 
grations of this form result in polylogarithms, of which the two we require 
are, 



t 



LUz) = / dt 



Li 2 (t) 



t 



(179) 



Both are real for z < 1. Their expansions for small z and for large — z are, 



Li a («)— ^z + OCa 2 ) , Li 2 (^)^-iln 2 (-^)-^ + 0(-) 



180) 



Li 3 (z)^z + 0(z 2 ) , Li 3 (z) 



- ln 3 (-^) - — ln(-z) + O(-) .(181) 



6 



6 



Care must be taken to arrange things so that the arguments of Li 2 (z) and 
Lis(z) lie in the range z < 1 for which the function is real. This is simple to 
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accomplish for the dilogarithm Li 2 (z) but it requires some effort for Li 3 (z). 
Li 3 (z) derives exclusively from Z 2 , when changing to dimensionless variables 
converts the rj' integration to either / dt/txhi (1—t) or / dt/t xln(t) ln(l— t). 
The final integral only makes sense for < t < 1 so it involves no choices, 



dt ln W ln(l-t) =Li3(t) _ ln(t)L . 2(t) 



(182) 



If the dimensionless parameter t lies in the range < t < 1 we write the first 
integral as, 



dt 



ln 2 (l-t) 
t 



-2Li 3 (l-t) + 21n(l-i)Li 2 (l-i) + ln 2 (l-£)ln(i) . (183) 



However, if t < we must use Landen's identity to re-express it as, 



dt 



ln 2 (l-t) 



-2Li3(t)-2Li3( T ^)-4Li 3 ( T i i ) + 21n( T ^)Li 2 ( T A) 

+41 <^) Ll2 (T^) +21n2 (T^) ln (l^) -( 184 ) 



It can happen that the appropriate choice to make between (11831) and (11841) 
depends upon whether x or r = ax is held fixed for large a, in which case we 
report the choice appropriate for holding x fixed at large a. 
Our results for the seven integrals are, 



-2xjf K = 2Li 2 (l-aiJx)-2Li 2 ( —— )+2Li 

r / 1+aHx 



1 ( 1 

-(l + — Hx 

2 v a 



2Li 2 [-(l + - + ^x)]+2Li 2 



i+l + aHx 



2Li, 



/ 1—aHx 

Wi- 



+ In 2 ( - ) - In (l + aHx) In [2 (- + Hx) ] + In 2 [2 ( 1 + - + Ex 

a 2(1- 



i+l — aHx' 



Hx 



-2xJ, 



GK 



2Li 3 ( 



—a 



+21n {-+Hx 



aHx+1 
Li 2 (-1)-Li 2 



■2Li, 



-2xJ, 



GK 
34 



A+- + 2Hx 

a 



-Li 



aHx 
aHx 
aHx—1 
1 — aHx 



T )+2Li 3 ( 



aHx 
aHx — 1 



-2Li, 



,2/1 tt \ , A — aHx 
+ ln 2 (- + Hx) ln(- — 

v a ) \l+aHx 

r-(l\ t-( l+ a Hx 
l2 \2J~ ^a+l + aHx 



(185) 



, (186) 
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1 , 9 /l+aHx\ 1 



1 
2 



1 + aHx 



+Li 2 — +- In 2 

^a+l-aHx J 2 v 

_i ln2(l+ I +fe) j + i ln(a+1 _ a&)+2(l+ I + ^ )lK _i+^|_ ) 

- 2 ( i 4-h in (d^)- 2 C+^) h^*) in ( i+a *) -( i87 ) 

a+l + aHx- 



H hifz 

a 



Ll *(^TT^)-Wl)+ln(. 



2 + 2aif x 



) 



+ \n(l + ^+Hx) ln( 



-Liofl) + ln 



a+l+aifa; 
2 + 2a#a; 
a+l—aHx 



1 



-if a; 



Ll2 C+l-aiix) 



a 



2+2aHx 
a+l—aHx 



+\n(l---Hx) ln(l+aHx 



1 \, / 1 \, /a+l—aHx 
1 + Hxj \n(l + — Hx) In 



a 



(188) 



2x^J — 



H hiix 

a 



Ll2 (a+l + afi:J U ^a 



2 + 2aHx x 
+ l + afW 



i+l + aifx' 

, /a+l + aHx\ , / 1 rr \ n / rr \1 _ _ . /a+l + aiirrx 
- ln (^T2^r) + ' n ( 1+ a + ^)' n ( 1+ai/:C )r 2L,2 (^^) 

(7 + 1- (lH.i\ 1 ,, . , /., ! ,. Nl I ' //.T\j 



-2Li 2 ( 



2a 



+ (l---ifx) jln 



1 1 

+ (l--+ii"x) ln(l--+ifa;) ln( 

\n(l---Hx) In 



a+1— aHx 
2+2aHx 



a+l+aHx 
2+2aHx 
a+l — aHx- 



-2xJ 9 GK 



2, /a+1 — aHx\ ^, . . , /a+i — a/ircx 
/-, x f T • fl-aHx\ /1\ / 1 + aHx \ 

-Uj-^+^^)-\n(2)\v(l + aHx)+\\n 2 (l+aHx)\ , (190) 



(191) 



It remains just to give the expansions for large a at fixed x, and at fixed 
r = ax. In both cases the leading term goes like ln(a). For fixed x we find, 



- 2xJ[ 



GK 



-4\n(2Hx) ln(a) +0(1) , 



(192) 
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-2xJ 2 GK 

— 2x l J^ 56 
2x J7^g 

-2xJ 9 GK 
-2xJ GK 



x ln(a) + 0(1) , 

41n(#x)ln(a) + 0(1) , 
2Hx-2(l-Hx) \n(l-Hx)] ln(a) + 0(1) 
2Hx + 2(l + Hx) ln(l+F:r)] ln(a) + 0(1) 

x ln(a) + 0(1) , 
-4 7 ln(a) + 0(1) . 



At fixed r the results are, 

-2xJ l GK = 

— 2xJ G i = 

—2xJ' 56 = 

2x^Ffg — 

-2xJ 9 GK = 
—2xJ 1Q 



x ln(a) + 0(1) , 

7T 2 /Fr-1 
— +2Li 2 ( 

-41n(l+ffr) ln(a) - 
x ln(a) + 0(1) , 
x ln(a) + 0(1) , 
x ln(a) + 0(1) , 
x ln(a) + 0(1) . 



ln(a) + 0(1) 
0(1), 



(193) 
(194) 
(195) 

(196) 

(197) 
(198) 



(199) 
(200) 

(201) 
(202) 
(203) 
(204) 
(205) 
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